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1.  STATEMENT  OF  PROBLEM  STUDIED 


Observed  and  modeled  concentration  fluctuations  in  smoke  plumes 
in  the  atmosphere  are  quite  large,  with  the  standard  deviation 
at  least  as  large  as  the  mean  "c".  There  is  interest  in  concentration 
fluctuations  because  of  their  importance  in  assessing  environmental 
and  toxicological  effects,  in  evaluating  models  for  predicting  mean 
concentrations,  and  in  determining  the  response  of  remote  sensors 
pointed  towards  smoke  plumes.  During  the  past  ten  years  there  has 
been  a  large  increase  in  research  on  this  subject  because  of  the 
concerns  listed  above,  because  of  advsn  es  in  the  development  of 
instruments  than  can  measure  short-term  concentration  fluctuations, 
and  became  of  increases  in  speed  and  storage  capabilities  of  computer 
that  some  researchers  are  using  to  run  models  of  concentration 
fluctuations. 

The  purpose  of  the  research  performed  under  this  contract  is  to 
review  existing  models  and  data  sets  on  concentration  fluctuations, 
and  use  this  information  to  develop  and  test  a  simplified  model. 
Emphasis  is  on  the  use  of  U.S.  Army  data  to  evaluate  the  model.  As  a 
result  of  the  model  evaluation,  research  required  for  model 
improvement  is  recommended  and  steps  taken  to  begin  some  of  this  new 
research. 

2.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

This  research  project  can  be  divided  into  three  tasks: 

e  Review  of  existing  models  and  data, 

e  Development  of  simple  model  and  evaluation  with  U.S.  Amy 
data. 

e  Development  of  strand  model. 

The  first  two  tasks  were  completed  during  the  first  20  months  of  the 
project  and  most  of  the  results  are  reported  in  the  publications  that 
are  listed  in  Section  3.  Abstracts  of  the  journal  articles  are  given 
here,  and  the  reader  is  referred  to  the  original  articles  for  details 


The  last  task  has  been  underway  since  September,  1984,  and  is  not 
completed.  It  is  briefly  reviewed  below  and  a  more  detailed  progress 
report  is  given  as  Appendix  A. 

2.1  Review  of  Existing  Models  and  Data 


CONCENTRATION  FLUCTUATIONS  IN  A  SMOKE  PLUME 
Steven  R.  Hanna 

Atmospheric  Environment.  18  (6) .  1091-1106 
(First  received  5  September  1983  and  received 
for  publication  17  January  1984) 

Abstract  -  Previous  research  results  are  reviewed  and  used  to  derive  a 
new  set  of  analytical  formulas  for  predicting  concentration 
fluctuations  in  smoke  plumes.  The  meandering  plume  approach  and  the 
internal  plume  approach  are  compared.  Some  simple  models  are  tested 
with  field  and  laboratory  data  sets,  showing  that  several  aspects  of 
the  data  (e.g.  the  concentration  fluctuations  on  the  plume  axis)  are 
reasonably  well  simulated.  However  there  is  much  room  for 
improvement,  since  the  models  have  some  fundamental  disagreements  and 
must  more  testing  with  data  should  take  place.  In  particular,  a  good 
knowledge  of  the  Lagrangian  time  scale  is  essential  for  predicting 
concentration  fluctuations. 

2.2  Development  of  Simple  Model  and  Evaluation  with  U.S.  Army  Data 

THE  EXPONENTIAL  PROBABILITY  DENSITY  FUNCTION  AND 
CONCENTRATION  FLUCTUATIONS  IN  SMOKE  PLUMES 
Steven  R.  Hanna 

Boundary  Layer  Meteorology.  29 ■  361-375. 

(Received  in  final  form  31  May,  1984) 

Abstract  -  Observations  of  1-s  average  concentration  fluctuations 
during  two  trials  of  a  U.S.  Army  diffusion  experiment  are  presented  and 
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compared  with  model  predictions  based  on  an  exponential  probability 
density  function  (pdf).  The  source  is  near  the  surface  and 
concentration  monitors  are  on  lines  about  30  to  100m  downwind  of  the 
source.  The  observed  ratio  of  the  standard  deviation  to  the  mean  of 
the  concentration  fluctuations  is  about  1.3  on  the  mean  plume  axis  and 
4  to  5  on  the  mean  plume  edges.  Plume  intermittency  (fraction  of 
non-zero  readings)  is  about  50%  on  the  mean  plume  axis  and  10%  on  the 
mean  plume  edges.  A  meandering  plume  model  is  combined  with  an 
exponential  pdf  assumption  to  produce  predictions  of  the  intermittency 
and  the  standard  deviation  of  the  concentration  fluctuations  that  are 
within  20%  of  the  observations. 

In  addition  to  the  above  paper  that  was  published  in  Boundary 
Layer  Meteorology,  three  other  papers  have  been  prepared  in  which  the 
simple  model  is  compared  with  observations,  two  of  these  can  be  found 
in  conference  proceedings  (CROC  Conference,  AMS  Conference),  and  the 
third  is  to  be  published  soon  in  Atmospheric  Environment.  The  first 
two  present  additional  statistics  from  the  U.S.  Army  Smoke  Week  III 
experiment,  including  observed  probability  density  functions  and 
moments  of  the  distributions  of  instantaneous  concentration,  plume 
width,  plume  centroid,  and  cross-wind  integrated  concentration.  The 
latter  paper  compares  the  new  model  with  observations  of  ground  level 
concentrations  resulting  from  tracer  releases  from  a  tall  stack  in 
Deardorff  and  Willis*  convective  tank. 

2.3  Development  of  Strand  Model 

The  strand  model  is  a  new  concept  that  is  under  development  and 
has  not  yet  been  published  in  a  journal  or  report.  A  preliminary 
description  of  the  model  is  in  Appendix  A.  The  purpose  of  this 
research  is  to  explicitly  account  for  concentration  fluctuations  at 
very  small  scales.  This  model  assumes  that  the  source  emits  strands 
of  polluted  material,  which  diffuse  due  to  molecular  motions  and  split 
when  their  diameters  double.  If  polluted  strands  are  close  enough 
together,  they  may  entrain  each  other  when  they  split.  The  whole 
packet  of  strands  also  diffuses  at  a  rate  comparable  to  that  of 


an  instantaneous  Gaussian  puff.  When  combined  with  a  model  for  the 
meander  of  the  instantaneous  plume,  the  model  gives  predictions  in 
agreement  with  wind  tunnel  data  of  Fackrell  and  Robins.  The  model 
produces  a  probability  density  function  of  the  concentration 
fluctuations  at  a  point,  and  can  account  for  the  effects  of  source  and 
sampler  size  and  sampling  and  averaging  time. 
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4.  LIST  OF  MEETINGS  ATTENDED 


Feb.  20-21,  1983,  White  Sands  Missile  Range. 

Concentration  fluctuation  data  sets  were  discussed  with  W. 
Ohmstede,  R.  Sutherland,  and  others  at  the  Atmospheric  Science 
Laboratory.  A  seminar  on  models  of  concentration  fluctuations 
was  given.  Plans  were  made  for  a  possible  smoke  experiment  to 
take  place  this  fall. 

June  15-17,  1983,  ARO  Quail  Roost  Workshop  on  Aerosol  Dispersion  in  the 
Atmospheric  Surface  Layer. 

At  this  workshop  current  U.S.  Army  research  projects  were 
discussed  and  recommendations  made  for  future  research  priorities. 

November  17,  1983,  Aberdeen  Proving  Ground,  MD. 

A  paper  entitled  "Concentration  fluctuations  in  smoke  plumes  from 
near  surface  releases"  was  presented  at  the  1983  Scientific 
Conference  on  Chemical  Defense  Research.  This  research  project 
was  discussed  with  D.  Sloop  and  R.  Saucier  of  CSL. 

December  5  &  6,  1983  and  February  3,  1984,  Boulder,  CO. 

Dr.  Hanna  attended  a  meeting  of  the  Large  Eddy  Simulation  Working 
Group,  sponsored  by  ARO.  Research  plans  for  this  project  were 
discussed  with  W.  Bach  of  ARO  and  W.  Ohmstede  of  ASL. 

June  25 ,  1984 ,  Aberdeen  Proving  Ground ,  MD . 

Dr.  Hanna  attended  the  1984  CRDC  Scientific  Conference  on 
Obscuration  and  Aerosol  Research  at  Aberdeen,  Proving  Ground,  MD, 
where  he  presented  a  paper  entitled  "Characteristics  of  Observed 
Concentration  Fluctuations  during  Smoke  Week  III. 


October  19,  1984,  Portland,  OR 


Dr.  Hanna  presented  a  paper  entitled  "Observed  and  Modeled 
Concentration  Fluctuations  in  a  Small  Smoke  Plume"  at  the 
AMS/APCA  Fourth  Joint  Conference  on  Applications  of  Air  Pollution 
Meteorology. 

October  23-25,  1984,  Kiawah  Island,  SC. 

Dr.  Hanna  attended  the  DOE/AMS  Workshop  on  Model  Evaluation  in 
Kiawah  Island,  SC,  where  he  presented  a  talk  on  the  results  of 
this  project  as  they  affect  model  evaluation,  and  discussed 
progress  with  Dr.  Bach  of  ARO. 

November  7,  1984,  Raleigh,  NC 

Dr.  Hanna  visited  Dr.  Bach  and  Dr.  Flood  at  ARO  in  RTP ,  NC. 

December  5,  1984,  White  Sands  Missile  Range,  NM. 

Dr.  Hanna  visited  ASL  at  White  Sands  Missile  Range,  NM,  where  he 
presented  a  seminar  on  his  research  and  discussed  this  project 
with  R.  Meyers  and  other  ASL  scientists. 

5.  LIST  OF  PARTICIPATING  SCIENTIFIC  PERSONNEL 

Dr.  Steven  R.  Hanna 
Dr.  Robert  J.  Yamartino 
Mr.  Jonathan  Pleim 
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1.  Introduction 


The  notion  that  pollutants  may  be  emitted,  dispersed  and 
described  mathematically  as  though  the  material  were  confined  to  small 
packets  of  size  X  has  been  described  by  Csanady  (1973) ,  Chatwin  and 
Sullivan  (1979),  Venkatram  (1983),  and  recently  by  Chatwin  (1984). 

The  prime  motivation  for  such  a  formulation  is  that  higher  moments  of 
the  concentration  distribution,  including  the  variance,  can  be  readily 
computed  since  the  problem  reduces  to  that  of  "counting"  statistics. 

The  size  of  the  packets,  X,  is  generally  chosen  to  correspond 
to  the  conduction  length  scale  where  molecular  diffusion  is  the  only 
operative  dilution  mechanism;  however,  most  treatments  then  either 
ignore  the  effects  of  molecular  diffusion  or  assume  that  packet  growth 
by  molecular  diffusion  results  only  in  the  entrainment  of  clean  air. 

In  this  paper  we  assume  that  pollutants  emerge  from  a  source  in 
long  thin  strands,  rather  than  small  three-dimensional  packets,  and 
that  these  strands  subsequently  grow  by  diffusive  processes  until  a 
time  td>  whereupon  they  split  into  two  strands  that  then  repeat  the 
above  mentioned  process.  Consideration  of  strand  entrainment  of  both 
clean  and  polluted  air  leads  to  a  differential  equation  for  the  time 
development  of  the  spatial  density  of  polluted  strands.  The  solution 
of  the  resulting  first-order,  non-linear,  differential  equation 
combined  with  counting  statistics  then  permits  a  determination  of 
ensemble  average  concentration  statistics  as  would  be  measured  by  a 
finite  aperture  detector  possessing  a  time  response  that  is  extremely 
rapid  compared  with  all  relevant  time  scales. 

Model  predictions  for  plume  intermittency  are  then  compared  with 
the  wind  tunnel  observations  of  Fackrell  and  Robins  (1982) .  Further 
assumptions  about  the  statistical  nature  of  the  strands  then  permits 
total  and  conditional  (i.e.,  C>0)  concentration  statistics  to  be 
evaluated.  Additions  to  the  theory  to  include  plume  meander,  time 
averaging,  and  finite  response  detectors  are  discussed. 


A -I 


2.  Strand  Theory 


2.1  Development  of  the  Governing  Equations 

Rather  than  postulating  a  non-linear  differential  equation  for 

the  concentration  variance,  as  is  done  by  Lewellen  and  Sykes  (1983), 

we  begin  with  the  simple  notion  of  pollutants  confined  to  thin  strands 

of  diameter  \  that  are  dispensed  by  turbulent  eddies  and  also  grow 

in  size  until  shear  forces  tear  the  strand  in  two.  Such  an  approach 

leads  naturally  to  the  computation  of  higher  moments  of  the  observed 

concentration  distribution  including  the  variance.  The  actual  size  of 

these  strands  will  correspond  to  the  conduction  length  scale,  which 

_3 

Chatwin  and  Sullivan  (1979)  estimate  to  be  of  order  10  meter; 

however,  specific  choice  of  a  numerical  value  for  X  is  not  a 

prerequisite  for  development  of  the  formalism. 

Consider  a  source  of  diameter  dg  emitting  at  rate  Q  (mass/time) 

into  a  flow  field  of  velocity  u.  Assuming  that  the  emissions  are  well 

mixed  across  the  source  aperture,  one  might  imagine  (Figure  1)  that 

the  pollutant  emerges  from  the  source  in  thin  strands  of  diameter  X 

and  containing  a  concentration  C  .  In  such  a  case  the  number  of 

o 

strands,  Ns,  emerging  would  be 

N  =  (d  /X)2  (1) 

s  s 


and  the  concentration  in  the  strands  would  be 


C  =  Q/(7  d  2  u)  (2a) 

O  4S 

or 


C  =  Q/[2ir  a  (0)  a  (0)  u]  (2b) 

o  y  z 

in  the  more  familiar  notation  of  plume  modeling.  Instantaneous 
concentrations  C(x,y,z)  downwind  could  then  be  expressed  as 

C(x,y,y)  =  C  D(x,y,z) 
o 


A-2 


(3) 


which 


where  the  instantaneous  plume  dilution  function  D(x,y,z)  is  defined 
such  that  at  the  source  D(x,y,z)  =  D(0,0,0)  =  1.  The  dilution 
function  is  assumed  to  be  known,  through  means  of  a  numerical  or 
analytic  model,  but  is  given  for  the  Gaussian  plume  model  as 

a  (0)  a  (0) 

D(x,y,z)  =  gy~X)~g~('X)  exP  2  (y/oy(x)>2  ~  2  (z/oz(x))2l  •  (4) 

y  z 


where  a  (x)  and  a  (x)  are  the  horizontal  and  vertical 
y  z 

instantaneous  plume  width  standard  deviations  respectively  as  a 
function  of  downwind  distance  x.  Substituting  Eqs.  (2b)  and  (4)  into 
(3)  gives  us  back  a  Gaussian  plume  expression  and  indicates  that 
nothing  unusual  results  from  this  re-grouping  of  terms. 

Ignoring  the  minor  "defect”  introduced  by  the  assumed  Gaussian 

shape  at  the  source,  earlier  assumed  to  be  uniform,  we  note  that  when 

the  dilution  function  takes  on  a  maximal  value  of  unity  all  available 

strands  are  populated  with  pollutant,  as  is  consistent  with  our 

definition  of  the  number  of  polluted  strands  N  by  Eq.  (1). 

s 

Temporarily  ignoring  diffusive  breakup  or  splitting  of  the  strands, 
these  entities  would  simply  disperse  relative  to  each  other,  so  that 
at  some  downwind  receptor  sampling  strands  of  air,  a  fraction  D(x,y,z) 
would  contain  pollutant  and  a  fraction  1-D  would  contain  only  "clean 
ambient  air".  Thus,  we  may  think  of  D(x,y,z)  as  either  the  usual 
diffusion  function  or  as  the  polluted  strand  density  function 
representing  that  fraction  of  locally  available  strands  that  are 
polluted.  With  this  information  in  hand,  one  could  immediately  write 
down  expressions  for  mean  square  concentration,  concentration 
variance,  and  other  moments  of  the  distribution;  however,  the  results 
would  be  in  poor  agreement  with  experiments  as  the  all-important 
diffusive  growth  of  the  strands  and  subsequent  strand  breakup  has  been 
ignored . 


If  the  strands  experience  diffusive  forces  specified  by  the 

molecular  diffusivity  v,  they  will  double  in  areal  size  in  a  time, 

t . ,  given  as 
a 

t.  =  \2/(4irv)  (5a) 

a 

whereupon  the  strand  is  sheared  apart  by  turbulent  eddies  into  two 

strands  of  diameter  X.  Given  molecular  dif fusivities  of  order 

-5  2  -3  -2  -3 

10  m  /sec  and  X  =  10  m,  one  computes  t.  in  the  10  -  10 

d 

second  range;  however,  as  with  X,  the  precise  numerical  value  of 

t^  is  not  critical  at  this  point  except  to  note  that  it  will  be 

sub-second  and  thus  relatively  short  in  comparison  to  many  time  scales 

relevant  to  atmospheric  dispersion.  More  generally  we  may  merely  say 

that  growth  in  the  number  of  strands  is  associated  with  a  time  scale 

t  without  relating  t  to  t.  via  a  rigorous  relation  such  as 

a 

t  =  t  /ln2  ,  (5b) 

a 

that  would  result  if  we  took  the  doubling  concept  literally.  Allowing 

such  a  growth  in  the  number  of  polluted  strands  creates  an  interesting 

dilemma  that  will  prove  quite  useful  in  developing  a  differential 

equation  for  the  number  density  of  polluted  strands.  Imagine  that  a 

long  pipe  with  diameter  dg  were  connected  to  the  source  described 

earlier  and  that  the  flow  speed  u  was  such  that  the  travel  time  in  the 

pipe  was  much  greater  than  t.  N  strands  of  concentration  C 

s  o 

would  go  into  the  pipe  and,  since  the  pipe  is  completely  filled  with 

strands,  N  strands  of  concentration  C  would  come  out  the  far 
s  o 

end.  What  if  anything,  happened  during  the  journey  along  the  pipe? 
Diffusion  and  strand  splitting  did  in  fact  occur,  but  since  the  pipe 
was  already  as  full  as  it  could  possibly  be  with  a  polluted  strand 
density  of  D  =  1,  the  growing  strands  had  no  other  choice  but  to 
devour  existing  strands  with  the  result  that  no  net  growth,  g, 
occurred  in  the  polluted  strand  density.  Can  we  then  write  down  a 
plausible  differential  equation  describing  the  time  evolution  of  the 
effective  polluted  strand  density,  p(t),  defined  as 

p(t)  =  g(t)  D(t)  ,  (6) 
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where  the  dilution  function  D(x,y,z)  has  been  re-expressed  as 
D(t)  =  D(ut,y,z)  for  compactness?  Knowing  also  that  the  number  of 
strands  will  grow  freely  as  g  ~  exp(t/r)  if  the  polluted  strand 
density  is  low  (i.e.,  p  <<1)  we  may  write  the  differential  equation 

p’  3  p(l-p)/T  +  g  D'  ,  (7a) 

where  the  prime  indicates  the  time  derivative  operator,  as  the 
simplest  form  that  possesses  the  three  necessary  properties: 

a)  exponential  growth  of  polluted  strand  fraction,  p,  with 
time  scale  t  for  p  small 

b)  termination  of  growth  as  p  approaches  unity  via  the  term 
(1-p),  and 

c)  forced  dilution  created  by  changes  in  D(t)  in  time. 

Applying  the  chain  rule  to  Eq.  (6)  and  dividing  through  by  D,_  Eq.  (7a) 
may  be  alternatively  expressed  as 

g'  =  g(l-p)/T  (7b) 

which,  because  the  factor  (1-p)  corrects  for  new  strands  consuming 
or  recombining  with  existing  strands,  elucidates  the  meaning  of  the 
term  "net  growth”  function  being  applied  to  g(t). 

Equation  (7a)  for  the  polluted  strand  density  or  fraction  is  the 
principal  result  given  to  us  by  this  simple,  physically-based  model  of 
strand  behavior,  and  this  equation  may  be  used  in  conjunction  with  any 
dispersion  model,  analytic  or  numerical,  giving  rise  to  a  diffusion 
function  D.  It  is  also  worth  noting  that  since  p(t)  is  closely 
related  to  the  intermittency  factor  y  discussed  extensively  by 
Wilson  et.  al.  (1984),  many  useful  expressions  could  be  written  down 
by  inspection:  however,  we  shall  first  explore  the  theory  further. 

2.2  Closed  Form  Solution 

Without  the  forced  dilution  function  D' ,  Equation  (7a)  is 
identical  to  that  used  by  population  biologists  to  describe  population 
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growth  in  a  world  with  a  fixed  and  finite  food  supply  and  is  easily 
solved  analytically  by  means  of  the  substitution  y  *  1/p.  (Kaplan, 
1958).  Using  this  hint  it  quickly  follows  that  a  solution4*  to 
Equation  (7b)  is 


g(t) 


-t/T 

e  +■ 


t 

dt’  e 


(t’-t) 


/TD(t*)}“1 


(8) 


where  t'  is  simply  a  dummy  variable  for  the  time  between  release 
(i.e.,  t  a  0)  and  present  time  t.  The  only  thing  disturbing  about 
Eq.  (8)  is  that  the  path  between  source  and  end  point  at  time  t  is  not 
specified.  Figure  1  shows  how  many  different  paths,  equivalent  to  the 
strand  trajectories,  could  end  up  at  a  specific  point  and  that  some 
weighted  consideration  of  all  possible  paths  (and  thus  strand  density 
time  developments)  should  be  involved  in  the  correct  solution  of 
Eq  (7b).  A  plausible  solution  might  consider  the  straight  line  path 
as  the  mean  solution  pathway  with  the  reciprocity  theorem  suggesting 
that  D(t')  in  Eq.  (8)  be  replaced  by 


D(t* ) 


2iro  (t-t'  )o  ( t 
y  z 


•t’)  JJ  dy’dz' D(ut’ 


y+y' 

m 


z  +z' 
m 


•  «xp{- 


1 

2 


(yr,)2 


1 

2 


(z  -z’ )2 

tit 

a  2(t-t’) 
z 


} 


(9) 


where  y  =  yt'/t  and  z  =  zt'/t  specify  the  mean  path  and  y’,z' 
tn  m 

are  dummy  integration  variables.  While  Eq.  (9)  may  represent  a 
significant  theoretical  point,  it  appears  to  be  of  secondary 
importance  in  most  applications.  In  addition,  we  note  that  Eq.  (9) 
can  be  cumbersome  and  numerically  problematic  to  evaluate  even  for  a 
single  mean  path;  thus,  we  turn  to  numerical  integration  of  the  more 
well-conditioned  Eq.  (7b). 


■►Non-linear  differential  equations  may  have  multiple  solutions. 
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2.3  Numerical  Solution 


Numerical  integration  of  Eq.  (7b),  subject  to  the  initial 
condition  g(0)=l,  must  recognize  that  the  potentially  explosive  growth 
(i.e.,  exponential)  in  the  solution  could  lead  to  rapid  divergence  of 
the  numerical  solution  from  the  true  solution.  For  this  reason  we 
adopt  a  second-order  accurate,  time-centered  (i.e.,  Crank-Nicolson) 
procedure  to  provide  a  reasonable  compromise  between  accuracy  and 
stability.  The  finite  difference  form  of  (7b)  then  becomes 


Vl  *  *n  ♦  S  "ntl  -  8n>a-2  Dn+l8n+l  '  2  W 


where  B  *  j  At/x  and  At  is  the  time  step.  The  time  advanced 
value,  g  ,  appears  on  both  sides  of  the  equation,  and  solution  of 
the  resulting  quadratic  equation  for  g^+1  yields  the  familiar 


(11 


where  a  =  ~  B  D  ,  , 

4.  n+ 1  _ 

and  c  -  -gn  [1  +  B(l-|  SnDn)]  . 

Equation  (11)  is  then  marched  in  time  using  a  time  step  chosen  as  the 
lesser  of  O.It  and  0.1  D/D'. 

Figure  2  shows  several  crosswind  profiles  of  the  resulting  p 
function  at  several  values  of  downwind  distance  x=ut  using  parameters 
relevant  to  the  wind  tunnel  studies  of  Fackrell  and  Robins  (1982). 

The  most  interesting  feature  of  these  curves  is  the  variation  in 
crosswind  profile  shape  as  one  moves  to  different  downwind  distances. 
Profiles  with  a  fairly  small  plume-axis  value  of  p,  denoted  pQ, 
look  Gaussian  in  shape  while  those  with  PQ-1  show  a  broad 
flat-top  region  with  a  precipitous  drop  at  large  y/c^.  As  it  will 
prove  important  in  the  inclusion  of  plume  meander  effects  to  know  this 


v  v '  ■ — mr  % 


- -  v  ■-  ^  v  ”  ^  ^7 1  ■.  '  •  '  ~  '■. 1  ^ v 


crosswind  profile  function,  the  crosswind  behavior  of  the  analytic 

solution,  Eq.  (8),  was  considered  assuming  that  the  on-axis  values 

were  available  from  the  numerical  integration  procedure  described 

above.  For  the  particular  trajectory  that  follows  the  path 

yaco,zaco,it  can  be  shown  that 
y  y’  z  z 

S(ut)  •  expt-«  {e2  +  c2)] 

p(ut,y,  z)/p  *  - r - *■ - - -  (12) 

1  +  exp[-±  (e„  +  c)]  [S(ut)-1] 
i  y  z 


where  pQ  =  p(ut,0,0) 

and  S(ut)  »  et/T/g(t)  =  D(ut,o,o)et/T/p  .  The  function 

o  ^ 

S(ut)  is  just  the  ratio  of  what  g  could  have  been  (i.e.,  e  )  if 

its  growth  were  unhindered  to  what  it  actually  is  on  the  plume 

centerline  and  hence  it  will  be  referred  to  as  the  "growth  suppression 

factor"  S.  Equation  (12)  provides  a  reasonable  description  of  the 

crosswind  profiles  shown  in  Figure  2  with  the  exception  that  falloff 

_2 

is  too  rapid  beyond  when  p/pQ  *  10  ,  but  this  may  not  be  a 

serious  drawback  since  this  plume  fringe  regime  should  not  dominate 
most  practical  problems. 

2.4  Statistical  Measures 


Once  the  mean  probability  p  is  known  at  the  receptor,  the 
binomial  distribution  (Ross,  1972)  gives  the  probability  P(j),  that  j 
of  the  detected  strands  contain  pollutant,  as 


N  ! 


-  nr  r  3 


0,1 . 


2 

where  Hr  *  (d^/X)  for  a  receptor  having  an  aperture  d^. 

Thus,  the  probability  distribution  for  detecting  j  polluted  strands  is 
identical  to  those  used  for  predicting  the  probabilities  in  die 
(p  »  1/6)  and  coin  (p  a  1/2)  toss  games. 
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Since,  there  is  always  the  finite  probability 

N 

P(0)  =  (1-p)  r 

that  none  of  the  strands  contain  pollutant,  the  "intermittency"  factor 

Y  =  1-P(0),  (15 

discussed  extensively  by  Wilson  et  al.  (1984),  will  be  less  than  one 
even  for  an  ideal,  zero  concentration  threshold  detector  and  in  the 
complete  absence  of  plume  meander,  as  presently  assumed  here. 

Figure  3  shows  the  intermittencies  predicted  by  Eq.  (15)  for 
plausible  strand  properties  of  t  =  0.2  sec  and  \  *  0.4  mm  and 
superimposed  onto  wind  tunnel  measurements  of  Fackrell  and  Robins 
(1982).  The  curves,  one  for  each  of  the  five  elevated  source 
diameters  considered,  demonstrate  that  the  strand  theory  approach  can 

i)  yield  predictions  that  exhibit  reasonable  qualitative 
behavior  in  terms  of  source  size  and  downwind  distance 
dependence  and 

ii)  span  the  range  of  observed  values. 

However,  as  important  ingredients  (e.g.,  plume  meander,  finite 
detector  thresholds  and  averaging  times,  strand  concentration 
distribution  assumptions)  are  still  missing  from  the  formulation, 
results  depicted  should  only  be  taken  as  grounds  for  "cautious 
optimism. " 

The  moment  generating  function  for  the  binomial  distribution 
gives  the  ensemble  average  number  of  polluted  strands  sampled  as 


the  ensemble  mean  square  number  of  polluted  strands  as 


.2 


N  <N  -l)p‘ 
r  r 


V 


(17) 


as  well  as  any  higher  order  measure  jn.  Thus,  the  ensemble  variance 
in  the  number  of  polluted  strands  is  just 


<J.2  =  j2  -  J  =  »rP(l-p)  .  (18) 

so  that,  even  if  at  each  point  all  strands  have  the  same  concentration, 
a  finite  ensemble  concentration  variance  would  be  observed. 

Hypothesizing  that  at  the  receptor  each  of  the  j  polluted  strands 
has  the  same  concentration,  Cg,  the  particular  sample  will  have, 
concentration 

C(j)  =  C  (j/H  )  (19) 

s  r 


and  the  ensemble  mean  concentration  will  be 


C  =  E  P(j)C(j)  =  C  (j/N  )  =  C  p  (20) 

.  s  r  s 

However,  for  this  result  to  be  consistent  with  Eq.  (3)  it  must  be  true 
that 


C  =  C  /g(t) 


(21) 


Similarly,  the  squared  concentration 


C2(j)  =  tCs<j/Nr)]2  =  j2[CQ/<gNr)]2 


results  in  the  ensemble  mean  square  concentration 


C2  =  I  P(j)C2(j)  =  [(1-1/M  )p2  +  p/N  HC  /g)2 
.  r  r  o 

J 


and  variance  of 


2  2~1  2 

<*c  =  C*  -  C  =  }-p<l-p><C  /g)* 

r 

The  total  concentrations  fluctuation  intensity,  i,  can  then  be 
expressed  as 

2  2  -2  1 

i  -  <*c'c  *  ±-u/P  -i> 

r 

Wilson  et  al.  (1984)  discuss  the  logic  and  convenience  of 
decomposing  the  fluctuation  intensity,  i,  into  two  parts:  the  part 
arising  from  detecting  positive  concentrations  and  the  part  due  to 
samples  where  the  detector  sees  no  material.  The  resulting 
decomposition  is 


.2 

l 


.2 

l 

P 


/y  +  (1-y)/Y 


.2 

l 


where 


and  the  subscript  p  refers  to  that  fraction  y  of  the  samples  where 
the  receptor  conditionally  sees  positive  (i.e.,  finite  non-zero) 
signals.  This  decomposition  into  conditional  quantities,  denoted  by 
the  subscript  p,  is  also  convenient  for  consideration  of  the  effects 
of  allowing  strands  to  have  a  distribution  of  concentrations. 

2.5  The  Exponential  Strand  Concentration  Distribution 

In  the  previous  section,  all  strands  at  a  given  point  in  space 
were  assumed  to  have  the  same  concentration  Cs;  however,  a  receptor 
at  some  downwind  location  will  detect  strands  that  have  experienced  a 
spectrum  of  histories  as  oppose*  to  the  simplistic  notion  that  each 
strand  will  have  suffered  g(t)-l  bifurcations.  In  developing  a  theory 
of  strand  "ageing”  as  a  sequence  of  bifurcations,  analogy  (Ross,  1972) 
with  other  decay  distribution  examples  (e.g.,  radioactive  decay, 
probabilistic  failure  analysis)  suggests  that  the  strands  will  have  an 
exponential  age  distribution 

,  1  -r/a 

P(r)  =  -  e  (28) 


where  r  is  the  age  parameter  and  a  represents  the  average  strand 
age.  As  increasing  strand  age  corresponds  to  a  greater  number  of 
bifurcations,  yielding  successively  lower  strand  concentrations,  it  is 
reasonable  to  conjecture  that  these  concentrations  also  will  be 
distributed  exponentially  as 

p  /  i 

C  (r)  =  C  e  (29) 

s  o 

where  6  is  a  constant  +  describing  the  effect  on  concentration  of 
a  single  bifurcation.  The  average  strand  concentration,  which  Eq. 

(21)  tells  us  must  be  CQ/g(t),  is  then 


+  1  r 

A  value  of  i  =  ln2  would  yield  C  (r)  =  C  (-)  implying 

S  0  2 

that  each  effective  bifurcation  reduced  the  strand  concentration  in 
half,  but  such  a  strong  assumption  is  not  necessary. 
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c  ® 

C  =  —  f  dr  exp[-r(l/a  +  1/6)] 
3  a  0 


=  C  /(I  +  a/6)  =  C  /g(t) 
o  o 


which  therefore  constrains  the  ratio  a/ 6  to  be 


(a/6)  =  g(t)  -  1 


and  leads  to  the  same  ensemble  mean  concentrations  of 


C  =  CQp/g 


obtained  assuming  uniform  concentration  strands.  Perhaps,  not 
surprisingly,  the  average  age  parameter  a  has  a  time  dependence 
proportional  to  g-1.  At  the  source  g=l  and  the  distribution  becomes 
an  infinite  spike  at  r=0  (i.e.,  the  delta  function  6(r)) 
corresponding  to  the  emission  of  unaged  strands  of  concentration 
Cq.  It  is  also  noteworthy  that,  even  at  great  distances  from  the 
source,  there  will  be  a  finite  probability  of  detecting  an  undiluted 
strand  with  the  original  concentration  Cq. 

Computation  of  higher  moments  of  the  observed  concentration  must 
recognize  that  for  any  one  sample  only  j  of  the  strands  are 
populated  with  pollutant.  Thus,  the  mean  square  concentration  for  j 
polluted  strands  involves  the  j-dimensional  integral 


C2(j) 


C2 

-fj  '  dri  l  d 

N  aJ 

r  o  o 


r  j  r  -r, /6  -r_/6 

J  dr^ [e  1  +  e  2 


-C./6.2 
+  e  j  ] 


-Cr,  +  r_  +  ...  +  r. )/a 
e  1  2  j 


C  .  .2  , 

-o  j - 3 -  +  - 3_=J - j 

ip  1  +  2  •/*  (1  +  a/4)2] 


and  the  ensemble  mean  square  concentration 


C2  =  I  P(j)C2(j) 

j 


involving  the  sum  over  j  from  zero  to  Nr  becomes 


—  C 


c  = 


M  p 

—  [ — - - 

n2  ll+2a/6 

r 


N  (N  -l)p  +Np-Np 
+  r  r _ 

(l+a/6)2 


2 

after  invoking  the  expressions  for  mean  j  and  j  given  by  Eqs.  (16) 
and  (17).  Cancelling  N^p  contributions  in  the  second  term  and 
replacing  at/6  with  g-1  then  yields  the  final  results 


_fi _ 

2g-l 


(Kr-l)p2 

2 

g 


and 
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a 


2 
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_£ _ 

2g-l 


-  (p/g)2] 


(3: 


for  the  ensemble  mean  square  concentration  and  concentration  variance 
respectively.  The  corresponding  conditional  estimates  can  then  be 
immediately  written  down  as  times  the  corresponding  ensemble 
mean  values  given  by  Eqs.  (32)  and  (34),  since  excluding  the  j  =  0 
terms  from  the  sum  requires  renormalization  of  the  probabilities  P(j) 
by  the  factor  =  (1-P(0))”^  to  ensure  that  the  conditional 
probabilities  Pc(j)  obey  the  constraint 


Z  P  (j)  *  1 


j  *  1,2,. .  .N 


«r 


’T'- 


where  P  (j)  =  P(j)/y.  The  equivalent  conditional  variance  is  then 
c 

just 


2 

a 

cp 


-  (p/g)2(H  (l/y-D+l)  ] 
r 


(36 


and  the  conditional  concentration  fluctuation  intensity,  i  ,  given  by 

P 


.2 

l 

P 


D(2g-1) 


-  N  (1/y-D-l ] 
r 


(37 


Far  from  the  source  g>>l  and  D<<1  so  that  the  first  term  dominates  and 
1/2 

i  *  [y/(2N  D) ]  ,  revealing  the  increasing  fluctuation 

P  r 

intensity  as  one  moves  away  from  the  plume  centerline. 

3.  Inclusion  of  Additional  Effects 


The  theory  development  to  this  point  has  focussed  on  the  problem 
of  determining  ensemble  statistical  measures.  Experimental  evaluation 
of  these  measures  would  require  repeated  sampling  at  fixed  points 
relative  to  the  centerline  of  the  instantaneous  plume.  While  perhaps 
not  an  impossible  experimental  feat,  it  is  of  greater  relevance  to 
adapt  the  theory  to  the  case  of  a  plume  meandering  back  and  forth 
across  a  receptor  fixed  in  space  and  accepting  samples  without  regard 
to  the  plume’s  position. 

In  the  case  of  the  Gaussian  plume  model  example  we  are  able  to 
build  upon  the  meandering  plume  formalism  developed  by  Gifford 
(1959).  Defining  a  coordinate  system  where  the  receptor  is  located  at 
a  point  (y,z)  relative  to  the  mean  centerline  position  of  the  plume, 
the  position  of  the  receptor  relative  to  the  plume  centerline  at  any 
instant  is  (y-y’.z-z’)  and  occurs  with  probability 


P(y’.z’)  = 


2iro 


exp [-  “  2(Z  /<Jifi)>2^ 
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neglecting  the  influence  of  ground  reflections.  Referring  back  to 
Eq.  (32)  and  recalling  that  D  =  p/g  we  may  immediately  recompute  the 
mean  concentration  as 


c  -  co  n 


dy'dz'P(y* ,z'  )D(x,y-y*  ,z-z' )  =  C^Cx.y.z) 


(3 


where  D  is  given  by  Eq.  (4)  and  D_  has  the  same  form  as  D  but  uses 

2  2  1/2  T 

the  total  sigma,  a  =  (a  +a  )  in  place  of  the 

T  M 


instantaneous  sigma  measures  a,  referred  to  as  a  and  a 

y  z 


throughout  this  paper. 

Computation  of  the  mean  square  concentration  involves  a  similar 
integral  averaging  of  Eq.  (34)  and  can  be  written 


dy'dz*P(y* ,z' ) +  1 


(40 


where 


D  =  D(x,y-y' ,z-z’ ) 
p  =  p(x, y-y’ , z— z  * ) 
g  =  p/D 


and  p(x, y-y',z-z')  can  be  approximated  with  the  help  of  Eq.  (12). 

Given  the  complexity  of  Eq.  (12),  evaluation  of  Eq.  (40)  will  either 
require  additional  approximations  or  numerical  integration. 

Inclusion  of  the  effects  of  a  finite  detector  concentration 
threshold,  CT>  is  another  complexity  that  may  require  numerical 
evaluation  in  most  cases.  Referring  back  to  Eq.  (33),  the  generalized 
form  of  the  nth  moment  may  be  expressed  as 


cn(j) 


,  •  «  «  -r  /i  -r_/4  -r./4 

(C^/N^)"*  ^  I  dr1  J  dr^...  J  dr.fe  +  e  +  . . .  +  e  J 

°  r  o  o  o  J 
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-(r, .r.)/aa/  -r,/6  -c./6  -c./i 

12  j  0(e  1  -t-e  2  +...+e  j  -C_N  /C  ) 

t  r  o 


where  the  theta  function  6(x)  is  defined  as 


f  0  for  x  <  0 
eu>  '  1  1  for  x  >  0 


Need  for  the  6-function  disappears  as  C^O  and  can  be  replaced 
by  altered  integration  limits  for  =  1  or  2,  but  we  have  not  yet 
developed  a  general,  analytic  solution  for  Eq.  (41);  however, 
integrals  of  this  type  are  relatively  common  in  statistical 
mechanics.  Once  accomplished,  one  could  compute  the  weighted  sum  over 
j  followed  by  the  integral  averaging  suggested  by  Eq.  (40). 

Another  important  effect  of  the  6  function  in  Eq.  (41)  is  that 
the  n  *  0  normalization  moment  is  less  than  unity  so  that  the 
probabilities  given  by  Eq.  (13)  are  modified  and  become 


P’(j)  =  P(j)  C  (j) 


This  has  the  immediate  effect  of  reducing  the  value  of  y  from  that 
given  by  Eq.  (15)  to 


y  =  1  -  P(0)  -  Z  P(j)[l-C  (j)] 
j=l 


thus,  indicating  the  significance  of  instrumental  "details"  on  the 
outcome  of  concentration  fluctuation  experiments. 
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Perhaps  the  last  refinement  that  would  be  considered  essential 

for  application  to  experimental  data  is  the  correction  for  averaging 

time.  For  very  small  averaging  times,  t  ,  much  smaller  than  the 

a 

Lagrangian  time  scale,  t  ,  associated  with  the  principal,  meander 

L 

driving  eddies,  it  is  thought  that  the  temporal  smearing  is  equivalent 

» 

to  an  increased  receptor  aperture  diameter  df  defined  as 


=  d  + 
r 


,  .1/2 
(a  a  )  t 
v  w  a 


(44) 


where  a  and  a  are  the  relevant  lateral  and  vertical 
v  w 

turbulent  velocities.  This  then  has  the  effect  of  increasing  the 
number  of  strands  received  by  the  receptor  to 

2  1/2  2 

N*  *  (d’/X)  =  N  CH-Cor  a  )  t  /d  ]  (45) 

r  r  r  v  w  a  r 

which  in  turn  has  the  effect  of  lowering  the  magnitude  of  the 

concentration  fluctuations  as  given  by  Eq.  (35)  for  example. 

Extension  of  the  theory  to  longer  time  averaging  periods  remains 

to  be  worked  out  in  detail,  but  may  require  further  assumptions  about 

the  form  of  the  auto-correlation  function  and  t  /t. . 

a  L 

Venkatram  (1979)  provides  results  of  assuming  an  exponential  form  for 
this  function  and  Hanna  (1984)  discusses  other  treatments  of  this 
problem. 

Finally,  it  should  be  noted  that  the  strand  theory  can  easily  be 
extended  to  multiple  species.  This  extension  could  be  an  important 
step  in  modeling  the  influence  of  sub-grid  scale  inhomogeneities  on 
processes  involving  non-linear  chemical  reactions. 
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